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We interpret the spectrum of meson states recently obtained in non-perturbative lattice QCD 
calculations in terms of constituent quark- ant iquark bound states and states, called 'hybrids', in 
which the qq pair is supplemented by an excitation of the gluonic field. We identify a lightest 
supermultiplet of hybrid mesons with J^*^ = (0, f,2) ^,1 built from a gluonic excitation of 
chromomagnetic character coupled to qq in an 5- wave. The next lightest hybrids are suggested 
to be quark orbital excitations with the same gluonic excitation, while the next distinct gluonic 
excitation is significantly heavier. Existing models of gluonic excitations are compared to these 
findings and possible phenomenological consequences explored. 



I. INTRODUCTION 

Models of hadrons as bound states of minimal num- 
bers of "constituent" or "dressed" quarks [H [2] can be 
motivated by the description of the flavor and J^^*-^^ sys- 
tematics of experimental meson and baryon spectra [3]. 
In particular, the absence of meson states with isospin 
greater than 1 or magnitude of strangeness greater than 
1 suggests a picture in which mesons are constituent 
quark-antiquark states where these constituent quarks 
have the same flavor quantum numbers as the quarks 
in the QCD Lagrangian. Higher Fock-state configura- 
tions, such as qqqq, would in general give rise to fla- 
vor conflgurations not observed experimentally. Further- 
more, all well-established states have J^*^ within the set 
0-+,0++,l--,l+-,l++,2-+,2++..., aU of which are 
accessible to a fermion-antifermion pair with orbital an- 
gular momentum, L. This suggests that the constituent 
quark degrees-of- freedom have spin-^ and the ordering 
of states in the spectrum implies increasing energy with 
increasing L. There is some theoretical support for the 
idea that constituent quarks with effective mass (at low 
energy scales) of several hundred MeV arise out of the al- 
most massless quarks of the QCD Lagrangian in the pro- 
cess of spontaneous chiral symmetry breaking (see e.g. 

mi)- 

Although the gross features of the experimental spec- 
trum strongly suggests this qq structure assignment for 
most states, such simplicity looks peculiar viewed from 
the perspective of QCD, in which quarks couple strongly 
to a self-interacting gluonic flcld. From the generic 
properties of QCD we might expect to have states in 
which the gluonic field itself is excited and carries J^'~^ 
quantum numbers. Absent any 'valence' quark content, 
such a state is called a 'glueball' and these basis states 
would be expected to appear mixed into the spectrum of 
isoscalar mesons [7] . The addition of a constituent quark- 
antiquark pair to an excited gluonic field gives us what 
we term a 'hybrid' meson, where these states have the 
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fiavor quantum numbers accessible to qq. They do not 

simply augment the regular qq spectrum however since 

p c 

the excited gluonic field could carry Jg" " quantum num- 
bers other than 0^"*". The gluonic quantum numbers can 
couple to the qq quantum numbers to give rise to so- 
called 'exotic' meson J^*^, literally those not accessible 
to a qq pair alone. Observation of a state with quantum 
numbers in this set, , 0"' , 1 '', 2^ ,3 ^ . . ., is con- 
sidered a smoking gun signature for states beyond the 
simple qq assignment. 

Now since the constituent qq picture was largely moti- 
vated by the absence of such exotic J^'-^ states in the 
observed spectrum, we must question why they have 
not been observed. Possibilities include a large energy 
scale associated with gluonic excitations, that places the 
states in a high mass region that has not been explored. 
Alternatively, the mass scale may be modest, but the 
states may have production/decay characteristics suffi- 
ciently different from qq mesons that experiments have 
thus far missed them [5J |S] . 

At present the experimental situation surrounding ex- 
otic J^'-^ mesons is confused, with 1 states being been 
claimed and disclaimed in several final states (see the re- 
view in :9J), while other exotic J^*^ are essentially unex- 
plored. There is hope that in the current decade we will 
see new data from GlueX, CLAS12, BESIII, Compass 
and PANDA that, through properly constrained analy- 
sis, will indicate definitively the presence or absence of 
exotic J^'-^ states in the light-meson spectrum. 

From the theoretical side, estimates of hybrid meson 
properties have traditionally followed from models which 
proceed from an assumed form for the gluonic excitation. 
There are a spread of such models and they make diver- 
gent predictions for the spectrum that, in the absence 
of clear experimental results, cannot be tested. Lattice 
QCD offers the possibility of computing meson properties 
directly from QCD, through numerical computation of 
the QCD path integral under the controlled approxima- 
tion of a finite, discretised, space-time grid. While lattice 
calculations have made impressive progress in precision 
computation of the masses of the lightest hadrons 10], 
it is only very recently that extraction of excited state 
properties has become feasible pTMl5) . 



2 



In this paper we show that the lattice QCD spec- 
troscopy data presented in p!2Hl4] can be interpreted 
in terms of a spectrum of constituent qq states supple- 
mented by hybrid mesons, with both exotic and non- 
exotic J^'~' . The lightest supcrmultiplct of hybrid mesons 
is clearly identified for the first time in a framework di- 
rectly connected to non-perturbative QCD. We compare 
the systematics of the extracted spectrum with expecta- 
tions of existing models of gluonic excitations and begin 
to develop a QCD-motivated phenomenology that may 
be compared to experimental data. 



II. GLUONIC EXCITATION MODELS 

A number of models proposing possible forms that the 
gluonic excitation within a hybrid meson might take have 
been explored [T6Vl21j . Notably, they make differing pre- 
dictions for the spectrum of hybrid mesons, so there is a 
hope of eventually distinguishing which, if any of them, 
resembles reality. 

In the flux-tube model [17], the gluonic field in a me- 
son is assumed to form itself into a tube between the 
quark and antiquark. Motivation comes from the strong 
coupling limit of QCD [22 , explicit computation of the 
action density between well-separated static color sources 
[531 HI] and the corresponding static potentials [551 US] ■ 
The model confines quarks and gives rise to linear Regge 
trajectories for conventional mesons. Hybrid mesons are 
described in terms of transverse oscillations of the tube, 
and the high degree of symmetry possessed by the system 
leads to a large degeneracy for the predicted lightest hy- 
brid supermultipleti; (0, 1, 2)"+, (0, 1, 2)+-, f— , 
The model can be extended by including flux-tube break- 
ing dynamics that give rise to meson decay |27| . 

Bag model approaches [HI [19] confine quarks inside 

a cavity. Gluonic excitations correspond to gluonic field 

modes allowed by the boundary conditions on the cavity 

waU and it is found that a "TE" mode with Jg"^" = 1+- 

is lightest, significantly lighter than a "TM" mode with 
p c 

Jg" " = 1 . This gives rise to a lightest hybrid super- 
multiplet (0, 1, 2) '", 1 if the qq is in an S-wave. Per- 
turbative elaborations of the model allow for non-exotic 
J-f^C" Jiybi-id states to mix with qq and for states to decay 
to pairs of lighter mesons. 

Another alternative description of gluonic excitations 
proposes that they correspond to the addition to the qq 
pair of one or more massive constituent gluons. By anal- 
ogy to constituent quarks being massive excitations of a 
non-trivial vacuum state filled with qq pairs, constituent 
gluons would be quasi-particle excitations of a non-trivial 
gluonic vacuum state. Particular implementations of this 
idea include the Coulomb-gauge many-body formalism 
in which the transverse quasi-gluons are excitations of 
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a variationally obtained approximate non-perturbative 
vacuum state [5ll20l[2T]. 

In the simplest constituent gluon picture one can add 
a Jg^'^" = 1 transverse quasi-gluon in a relative S- 
wave with respect to a qq pair. This gives rise to a su- 
pcrmultiplct (0, 1,2)++,!^ if the qq is in an S-wave. 
This exactly resembles a qq P-wave and would be very 
hard to pick out of a spectrum unless some distinguish- 
ing decay characteristics could be observed. With the qq 
pair in a P-wave however, one has a large supermulti- 
plet 0~-,(l— )3,(2~)2,3~,0-+,l-+,2-+, featuring 
negative parity exotic states. Notably there is a light 
exotic - a prediction essentially unique to this model 

[iai2o]. 

There are arguments based upon the existence of re- 
pulsive three-body (qqg) forces [21] and the behavior of 
adiabatic potentials in QCD at short distances [3B], that 
in fact having the quasi-gluon in a P-wave with respect 
to the qq, coupled to total gluon spin Jg" " = l+~, is 
energetically favorable. In this case one has a lightest hy- 
brid supermultiplet of (0,1,2) '",1 with the qq pair 
in an S'-wave and, somewhat heavier, a larger supermul- 
tiplet of O+-,(l+-)3,(2+-)2,3+-,(0, 1,2)++ where the 
qq are in a P-wave. 

Although they make very different predictions for the 
spectrum, without experimental evidence for a pattern 
of hybrid mesons it is difficult to favor any one of these 
models above the others. It is here that we will take ad- 
vantage of recent improvements in spectroscopic calcula- 
tions in lattice QCD to provide solid theoretical evidence 
for a spectrum of hybrid mesons with both exotic and 
non-exotic J^'-^ . 



III. LATTICE QCD SPECTROSCOPY 

Lattice QCD offers a means to obtain hadronic quanti- 
ties directly from QCD by computing appropriate gauge- 
invariant correlation functions constructed from quark 
and gluon fields. By considering it on a finite grid, the 
QCD path integral can be Monte-Carlo sampled and cor- 
relation functions determined with finite statistical preci- 
sion. In principal one can, by performing multiple calcu- 
lations, extrapolate in decreasing lattice spacing, and in- 
cresing finite volume of the box - in this sense the lattice 
approximation is systematically improvable. There are 
technical challenges, notably the poor scaling of the com- 
putational cost with decreasing quark mass, which means 
that most calculations currently have to be performed 
with unphysically heavy up and down quark masses. 

The spectrum of hadrons follows most directly from 
two-point correlation functions, 

C,j{t) = (0|O,(i)O,(0)|0) 

where O are gauge-invariant combinations of quark and 
gluon fields (expressed via the parallel transporters, or 
"links" of the lattice) that have the desired hadron quan- 
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turn numbers. Detailed discussion of the methodol- 
ogy required to extract an excited state spectrum with 
determined continuum J^*^ quantum numbers can be 
found in [13] - here we summarise just the salient points. 
The quantities extracted from the calculations are the 
state masses, mn, and vacuum-operator-state "overlaps", 
= (n|Oi|0), that appear in a spectral decomposition 
of a two-point correlation function, 



Q,(t) = (0|O,(t)Oj(0)|0)=^ 



where the sum is over all eigenstates, |n), of the QCD 
Hamiltonian in finite volume with the same quantum 
numbers as the operators Oij. There is a variational 
method of analysis [551122] that analyses a matrix of such 
correlators built using a basis of operators, returning best 
estimates for masses and overlaps. 

In jT3] , a large operator basis for mesons was built us- 
ing fermion bilinears projected onto zero meson momen- 
tum, featuring up to three gauge- covariant derivatives^ 
(d = 0,1,2,3), 



{Jrmr; JDmD\J, M) 



where 



Vo = l,m = Am) ■ D , 

and where the three-derivative construction can be found 
in [T3|. The use of a circular polarisation basis, e (to), for 
the derivatives enables the use of simple 5*0(3) Clebsch- 
Gordan coefhcients for the angular momentum construc- 
tions. 

The quark fields featuring in these expressions do not 
just exist on a single site on the lattice - they are smeared 
over space in a gauge-invariant way. The particular form 
of smearing used is "distillation" [30] , which proves to be 
a highly efficient way to compute a large number of corre- 
lators with operators that sample dominantly the lowest 
energy modes of QCD relevant to low-lying hadron states. 
Furthermore all gauge-fields entering into the operator 
construction are smeared in a gauge-covariant manner 
known as "stout-smearing" [3T] that preserves transfor- 
mational properties of the original fields. The consistent 
use of gauge-invariant constructions is relevant since the 
computation is not performed in any fixed gauge. 

The calculations reported in |13| were performed on 
dynamical anisotropic Clover lattices with three flavors 



of quark, the lightest two of which are mass degen- 
erate and the third is tuned to describe the strange 
quark. The lattice spacing in the spatial directions is 
fls ^ 0.12 fm, and this proves to be fine enough to see an 
effective restoration of rotational symmetry at the scale 
of hadrons. The temporal direction has a finer spacing 
corresponding to d^^ ^ 5.6 GeV that gives an excellent 
resolution of the time-dependence of hadron two-point 
correlators. With the strange quark mass held fixed, 
light quark masses corresponding to pions of mass be- 
tween 400 and 700 MeV were considered. Spatial vol- 
umes were ~ (2.0 fm)^, (2.5 fm)'^ and no significant 
trends in volume dependence were observed in the ex- 
tracted spectra. This fact was used to argue that these 
calculations are not resolving the expected physics of ex- 
cited states as resonances in meson-meson scattering. In 
a finite volume one would expect a discrete and volume- 
dependent spectrum of meson-meson energy levels that 
were not observed. The spectrum that was observed was 
thus interpreted as 'single-hadron' states whose resonant 
nature cannot be explored without adding to the calcula- 
tion operators resembling pairs of hadrons. In this paper 
we will explore the possible bound-state compositions of 
these single-hadron states. 

Quantities extracted from the lattice calculation are di- 
mensionless, dimensionful quantities being scaled by an 
appropriate power of the lattice spacing - for example 
a mass would be extracted as atm, where, if at is deter- 
mined, we can scale-set and obtain a mass in MeV. Scale- 
setting in a calculation with unphysically heavy quarks 
is an ambiguous procedure since there is no experimental 
quantity that we can compare with to set the scale. We 
choose a simple but mass-dependent procedure to set the 
scale - all mass-dimension quantities are scaled as follows. 
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where (atTOn) is the dimensionless mass of the f2 baryon 
computed on this lattice and to^^'''^'*' is the experimental 
mass. 

We will focus initially in this paper on the spectrum 
results obtained at the heaviest 'light' quark mass, where 
all three quarks are at the strange quark mass and the 
corresponding pion mass is ~ 700 MeV. We will see that 
the qualitative features observed in the spectrum will 
change relatively little as the pion mass is reduced to 
~ 400 MeV. The spectrum of isovector states extracted 
is shown in Figure [l] where we see the benefits of using a 
large basis in the number of excited states extracted in a 
given J^^ and the range of J^^ considered. 



^ The "forward-backward" gauge-covariant derivative Dp = 
— 2igAij, makes construction of operators of defi- 
nite charge-conjugation relatively simple. On a cubic lattice 
the derivatives are implemented via parallel-transported finite- 
differences featuring the SU{3) matrices living on the links. 



IV. INTERPRETATION OF THE MESON 
SPECTRUM 

Of immediate importance is the presence in Figure [T] 
of a spectrum of exotic J^'^ mesons: a lightest 1 ^, sig- 
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FIG. 1: Spectrum of isovector mesons extracted in [13] with a pion mass of ~ 700 MeV. Masses (in MeV) are averages over two 
computed volumes, ~ (2.0 fm)'^, (2.5 fm)^. Box size indicates the statistical uncertainty on the extracted mass. Ellipses indicate 
that there may be heavier states in that J^'^ but we did not determine them reliably in [T3] . States are color coded by assigned 
I/-wave supermultiplet as described in the text. States colored grey could not be assigned to an L-wave supermultiplet. 
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TABLE I: qq ^^+^Lj, J^'^ supermuhiplets. 



nificantly lighter than positive parity states, 0"' , (2^ Y . 
As well as more 1 ^ states, much higher in the spectrum 
are and 3 states. We will argue that these exotic 
states can be explained as being hybrid mesons. 

Under examination, the spectrum of non-exotic J^'^ 
mesons shown in Figure [l] displays J^*^ patterns con- 
sistent with qq^^^^Lj structure (as presented in Table 
[l| . These patterns are identified in Figure [ij where we 
show apparently complete S,P,D, F-wave supermulti- 
plets^ , and a partial G-wave supermultiplet indicated 



by the presence of 4 ,4 states. There appear to 
be higher mass recurrences of several of these supermul- 
tiplets, which in analysis to be described later we will 
suggest could be radial excitations of qq. 

Once these assignments based upon degeneracy pat- 
terns are made, we can see a near-degenerate set of states 
with J^'^ = '",1 ,2 remains unassigned (colored 
blue and indicated with a star (*) in Figure[I]). No super- 
multiplet in Table |T] contains such a set, and we notice 
that the mass-scale of these states is very close to the 
lightest exotic state, the 1 '". We hypothesize that these 
states could be members of the lightest hybrid meson 
supermultiplet . 

To explore the hypothesis that we have a supermulti- 
plet of hybrid mesons embedded within a spectrum of qq 
states, we take advantage of the "overlap" information 
provided by the two-point correlator analysis described 
in the previous section. In the same variational analysis 
that determines the mass spectrum, nin, we also deter- 
mine the matrix-elements, Z" — (^n\Oi\0'), that encode 
the degree to which operator d "overlaps with" state n. 
To the extent that we understand the structure of our 
operators, we can infer some information about the in- 
ternal structure of the extracted states from the relative 
size of these overlaps. 



^ The terminology "supermultiplet" is a historical one, where fla- 
vor multiplets of a given J^'-^ were first identified and later dif- 



ferent J^'-^ having a common proposed bound-state structure 
were related. 
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A. Operator interpretation 

As an example, consider the simplest 1 operator in- 
cluded in our basis, '07iV'i which as a shorthand notation 
we call "p" . If we argue that the smeared fermion fields 
are capable of producing a single constituent qq pair, 
then writing out the quark-spin and momentum struc- 
ture of the spinor contraction, Ua-{Pq)jiVg-{pq), we realise 
that this operator will dominantly overlap with qq ^Si 
states. Overlap with qq ^Di comes from lower compo- 
nents of the Dirac spinors and is suppressed by powers 
of This "sub-leading" overlap can be further sup- 
pressed by considering the upper-component projected 
operator ■!/57ji(l — 7o)'0, which we call "pnr"- 

We must emphasise here that since the fermion fields 
are smeared over space we are not measuring the decay 
constant of a vector meson with the "p" operator. Rather 
we are sampling a smearing function weighted integral 
of the state's wavefunction 32] . A bound-state model 
description of this can be found in Appendix A. 

To extend our interpretative framework consider the 
simplest operator featuring a gauge-covariant derivative, 
i/j D ill), which we call (clq x -Djlj^) ^ and which trans- 
forms like J^*^ — 1 . Initially let us simply neglect 
the gauge-field term in the gauge-covariant derivative in 
our interpretation and use the ordinary derivative act- 
ing on the quark field to introduce an extra factor of the 
quark momentum. The resulting spinor contraction has 
overlaps of comparable size for both qq ^Si_ and qq ^Di . 

At the two derivative level, we have an example of an 
operator which has leading overlap with qq ^Di and sup- 
pressed qq'^Si, namely (pnr x D^jL2) ■ Interpreting 
the ordinary derivative part of the operator as giving rise 
to a factor Y^"^ (d) we see where the Z?-wave structure en- 
ters. 

The same basic procedure of neglecting (for interpre- 
tation purposes) the gauge-field part of gauge-covariant 
derivatives can be followed for all our operators, yield- 
ing qq ^^~^^Lj constructions for non-exotic J^*^, except 
that once we reach two-derivative constructions we hit 

[21 

a non-trivial case. One possible combination is -D vli = 

(l, TOi; 1, 7712 1 1, m) Z? D rn.2 1 wherc, since the Clebsch- 
Gordan coefficient is antisymmetric under exchange of 
(mi, 77x2), in the limit of ignoring the gauge-field part of 
the derivatives we have zero. This operator construc- 
tion is only non-zero precisely because the gauge-field is 
present and in fact is proportional to chromomagnetic 
components of the field-strength tensor, F. Since one 
must have a non-trivial gauge-field structure to have a 
non-zero value of F we propose that overlap with this 
operator indicates hybrid-like character. 

As a concrete example of such an operator we have 
which is an operator with exotic J^'~^ — 
The structure of this operator would appear to be 
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that of a qq pair with 5=1 and L — Q (from pnr) along 
with a non-trivial gluonic field having chromomagnetic 
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FIG. 2: Spectrum of 1 mesons (as in Figure [T]) along 
with relative size of overlaps onto operators, pnr ~ '^Si, 
Pnr X -Djl2 ~ ^Di and ttnr x -Djl^ ~ ^hybj. Histograms 
normalised such that the state with largest overlap across the 
entire excited 1~~ spectrum has mean value of 1. The lighter 
colored area at the head of each bar represents the statistical 
uncertainty. 



PC ^ 

character ( Jg " ' — V , from D 
of hybrid meson construction proposed in the bag model 
PEl \W\ and in the "P-wave" quasi-gluon picture [H] . In 
the lattice calculation we find that the lightest extracted 
1 ^ state has large overlap onto this operator. 

Now let us return to the spectrum of non-exotic J^'-^ 
states. We can ask how well the assignments of L-wave 
supermultiplets made on the basis of mass degeneracy 
patterns match with the overlaps onto operators having 
dominant qq '^^^^Lj character? As an example, in Figure 
[2] we present the 1 spectrum and the overlaps onto 
three characteristic operators: pnr, (pnr x £'jj^2) 
and (ttnr x -Djii)''"^. These three operators are dom- 
inantly qq^Si, qq^Di and, according to the proposal 
above, hybrid with quarks in a spin-singlet. The ground 
state we observe to be dominantly qq ^Si , and similarly 
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FIG. 3: Z values in lattice units. Left pane: lightest P-wave supermultiplet, 1^^,(0,1,2)^"*" overlaps onto {7r,p}NR x -Djlj^. 
Middle pane: lightest D-wave supermultiplet, 2 '',(1,2,3)"" overlaps onto {tt, p}nh, x -Djl2. Right pane: lightest hybrid 
multiplet, 1 , (0, 1, 2) ^ overlaps onto {tt, p}nr x -Djii- 



the first excited state, which we interpret as a radial exci- 
tation, having the same angular and spin structure as the 
ground state^. The second excited state, which we earlier 
assigned to a D-wave supermultiplet indeed is dominated 
by qq ^Di. It is notable that there appears to be rela- 
tively little mixing between qq ^Si, qq^Di basis states, 
even though the resulting eigenstates are very close in 
mass. The third excited state is the 1 state that could 
not be placed into an L-wave supermultiplet on the ba- 
sis of degeneracy patterns. We find that it has dominant 
overlap with the operator (vtnr x -Djlj^) that we have 
characterised as being a spin-singlet hybrid meson. 

Very similar conclusions can be drawn in the 2 ^ chan- 
nel^, where the ground state is found to be dominantly 
qq ^I?2 7 while the first excited state overlaps strongly onto 

(pnr X i'jl]^) ^ and which we thus propose is a hybrid 
meson with quarks in a spin-triplet. 

The ^ channel is not quite so clearly delineated - 
the ground state has a large overlap onto operators of 
qq ^Sq character, but also onto a hybrid (pnr x D^]^^^ 
operator. That the pion might have behavior not en- 
tirely compatible with constituent qq structure is to be 
expected given its role as the pseudo Goldstone boson of 
spontaneous chiral symmetry breaking, and we will not 
explore it in detail here. The first excited state appears 
more conventional, with overlap onto qq ^Sq operators 
and suppressed overlap onto (pnr x The sec- 

ond excited state has the largest (pnr x over- 
lap, but unlike the 1 ,2 *" hybrid candidates, it also 
has significant overlap onto qq operators. This may indi- 
cate a larger degree of mixing with qq basis states than 
is present in the other non-exotic channels present in the 



* We should point out here that a radial excitation can have an 
'accidentally' small overlap with an operator whose angular de- 
pendence accurately characterises the state's structure - see Ap- 
pendix A for an explicit demonstration 

^ see Figure I4I 



supermultiplet. 

In identifying a supermultiplet, one expects the Z val- 
ues from a common operator projected across the allowed 
J^^ to be similar. For example, pnr x -Djl^ can be pro- 
jected into J^'" — (0, 1, 2)++ using appropriate Clebsch- 
Gordans and the corresponding spin-singlet comes from 
ttnr X i^jl^- The normalisation is such that one would 
expect the Z value extracted to be common to all mem- 
bers of the supermultiplet. We show the extracted Z val- 
ues for the lightest proposed qq P-wave supermultiplet in 
the left-hand pane of Figure [3] where we see reasonable 
agreement, supporting the assignment. In the middle 
pane of Figure [3] we show a similar pattern for the light- 
est proposed qq D-waye supermultiplet. The right-hand 
pane shows the pattern for our proposed lightest hybrid 

multiplet using the operators {{TT,p}m. x D^jLi) ^ ^ that 
strongly suggests these four states have a common struc- 
ture^. 

We caution the reader at this point that only qualita- 
tive judgements can be drawn from the value of the Zs 
since these matrix-elements determined with lattice regu- 
larisation require renormalisation to be compared to the 
continuum theory matrix-elements. That renormalisa- 
tion can mix operators, including those of differing mass 
dimension, although with a relatively fine lattice spac- 
ing, an improved action, and smeared fields we do not 
anticipate such mixing being large. 

We have now demonstrated the main result of this 
paper, the first identification of the lightest hybrid su- 
permultiplet within QCD. As well as a single exotic 
1 *" state, there are three non-exotic quantum numbered 
states, 1 , 2 that differ from the corresponding 
qq constructions by having the opposite quark-spin (sin- 
glet versus triplet). 

The reader may be concerned that the results pre- 



A model-dependent demonstration using a constituent gluon con- 
struction is given in Appendix B. 
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sented thus far are for a version of QCD in which the 
three 'hght' quark flavors are all at the strange quark 
mass, and this world may have little to do with the ver- 
sion of QCD manifested in this universe which has two 
flavors being much lighter than strange. The possibil- 
ity of significant changes under reduction of the light 
quark mass can be explored in lattice QCD calcula- 
tions with lighter quark masses and calculations down 
to rriTr ~ 400 MeV were reported in [T3] . The same anal- 
ysis in terms of qq L-wave supermultiplet degeneracies 
and operator overlaps has been carried out and the qual- 
itative features remain unchanged. Relevant to the par- 
ticularly important case of the lightest hybrid supermul- 
tiplet, we show in Figure |4] the quark mass dependence 
of the 1 ,0 ^ and 2 spectra and operator overlaps. 
We see that qualitatively there is relatively little change, 
except possibly in the degree of gg-hybrid mixing in the 
(0,2) channels. The mass splitting within the super- 
multiplet grows with decreasing quark mass as shown in 
Figure [5| but the lightest negative parity hybrid states 
remain systematically below the positive parity exotics 
and we argue that it remains sensible to describe them 
as a supermultiplet^. 

In the remainder of the paper we will be concerned 
with extending our description of the spectrum, identi- 
fying higher mass hybrid states, comparing with models 
of gluonic excitations and finally with developing a phe- 
nomenology of hybrid mesons that is directly motivated 
by results of QCD calculations. 



B. Heavier hybrid mesons 

Having identified what we propose is the lightest su- 
permultiplet of hybrid mesons, we now move on to con- 
sider the presence of heavier hybrid meson states in 
the extracted lattice QCD spectrum. The first obvi- 
ous candidates in the ~ 700 MeV calculation are 
the exotic 0"' and 2"' states above 2.4 GeV. The 
0^ state is found to have considerable overlap with 



the (ai X Djl]^)^ " operator, which has the same glu- 

[2] 

onic structure {Dj!^^) as the hybrids previously iden- 
tified, but which places the spin-triplet qq pair in a 
relative P-wave (ai ~ V'757iV' ~ IQ^Pi)- This sug- 
gests that there might be an excited hybrid supermul- 
tiplet with the same chromomagnetic gluonic excita- 
tion, but also excited in quark angular momentum to 
Lgg = 1. Such a supermultiplet would be expected to 
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^ the large mass of the 2~+ state may be due to significant mixing 
with a heavier qq ^D2 basis state as suggested by the overlaps in 
Figure [1] 
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FIG. 4: Mass spectrum and histograms as in Figure [2] as a 
function of quark mass expressed via the pion mass. 



spectrum with relatively little splitting between them. 
In order to find non-exotic members of this supermulti- 
plet embedded within the qq states one should examine 
the overlaps for a wider set of operators than we have 
so far considered including those constructed from three 
derivatives. We will not go into that level of detail in this 
paper. 

There remain further exotic states at higher mass. The 
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FIG. 5: Mass spectrum of lightest isovector hybrid supermultiplet along with lightest positive parity exotic states computed at 
four pion masses. 



heavier 1 state at about 2.6 GeV has a pattern of 
overlaps that is similar to the lightest 1 ^ and it thus 
may correspond to a 'radial excitation' in which the 
quark radial motion is excited on top of the same chro- 
momagnetic gluonic excitation. The heaviest 1 state 
shown in Figure [21 which has a significant overlap onto 

[21 

ttnr X D j^^, could conceivably be another member of 
a radially excited hybrid supermultiplet - the remaining 
(0,2) members are not identified in this calculation. 

The 3 ^ state above 3 GeV could be part of a super- 
multiplet having qq in D-wave coupled to the lowest chro- 
momagnetic gluonic excitation, or it could arise from a 
different gluonic excitation. One possibility would be S- 
wave qq with gluonic field transforming as Jg " ' = 2"' . 
At the level of detail we have in this calculation we can- 
not definitively state which is the case. 

There is an exotic state that appears at a very 

heavy mass above 3 GeV - such a state does not arise with 

a chromomagnetic gluonic excitation combined with qq 

and most likely signals the mass scale for a different type 

p c 

of gluonic excitation, possibly one having Jg" ^ — 1 . 

p c 

We note here that the energy ordering of Jg " " as 
1^ <1 <... is the same as is observed for gluelumps, 
the excitations of the gluonic field around a color-octet 
source in SU{i) Yang-Mills theory [26]. 

C. Isoscalar mesons and kaons 

So far we have identified J^'~' sztpermultiplets within 
the isovector meson spectrum but have not considered 
flavour partners in a flavour multiplet. We expect there 
to be kaonic and isoscalar states within flavor octets and 
singlets associated with the isovector states. Kaons have 
a denser spectrum described by the more limited 
quantum numbers, the flavored states not being eigen- 
states of charge-conjugation (or the G-parity extension). 
An immediate consequence of this is that there are no 
exotic kaons, and hence no smoking gun signature for 
hybrid kaons. Within simple qq models, the spectrum of 
kaons features states constructed from admixtures of op- 
posite C, for example in the axial (1"*") sector one would 
expect two low-lying states constructed from the basis 



states us ■^Pi(l+(+)), us -^Pi(I+(^'). Experimentally two 
states are found, the iiri(1270), i^i(1400), whose decay 
properties suggest that they are strong admixtures of the 
above basis states, with a mixing angle close to 45° [3]. 

Kaon spectra were extracted in the lattice calculations 
reported in (I3j using the derivative basis described ear- 
lier, but using both and J^^ operators together 
in the calculation. As such, by considering overlaps 
(n(J^)|O(J^(c^))|0) one can infer the degree of mixing 
between C = + and C = — basis states. In explicit com- 
putation we found this mixing to be very small, but this 
is likely to be a reflection of the unphysically heavy light 
quark masses used (the lightest is ~ 400 MeV) which 
do not give rise to as large a breaking of S'J7(3)f as is 
present in nature. In light of this, and the impossibility 
of exotic quantum numbered states, we will not discuss 
kaons any further. 

Isoscalar meson properties are somewhat more chal- 
lenging to compute in lattice QCD owing to the need 
to evaluate disconnected Wick contractions. Distillation 
[30] proves to be an efficient procedure to do this and 
in [14] we reported a calculation of the isoscalar meson 
spectrum at a single pion mass, TOtt ~ 400 MeV using the 
operator basis of [13] , doubled in size by including both 
light-quark (£ = ^ {uTu + dTd) ) and strange-quark bi- 

linears. The spectrum obtained in a (2.0 fm)'^ box is re- 
produced in Figure |6j The degree to which a given state 
is dominated by ^ or s can be estimated using the rela- 
tive size of overlaps onto operators. We parameterize the 
admixture by introducing a mixing angle, a, shown in 
Figure |6j assuming pairs of states (o, b) are orthogonal 
combinations of just two light-strange basis states: 

I a) = cos a I ^) — sin a I s) 
|b) = sin a I ^) + cosa|s). 

We note that the basic qualitative phenomenology ob- 
served in experiment is reproduced. In most J^*^ the 
flavor mixing is close to ideal, giving rise to almost pure 
[uu+dd}) and ss states. The exceptions to this are the 

^ channel where the rj, r( mixing being quite close to 
SU(?))p ideal octet-singlet is observed and the 1"'""'" chan- 
nel where phenomenological analysis of experimental ra- 
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FIG. 6: Isoscalar meson spectrum with ~ 400 MeV labeled by J^^ (reproduced from [T^]). The light-strange content 
of each state (cos^ a, siv? a) is given by the fraction of (black, green) and the mixing angle for identified pairs is also shown. 
Horizontal square braces with ellipses indicate that additional states were extracted in this J^*^ but were not robust. Grey 
boxes indicate the positions of isovector meson states extracted on the same lattice (taken from [13]). Pink boxes indicate the 
position of glueballs in the quark-less Yang-Mills theory [33j . The candidate states for the lightest hybrid meson supermultiplet 
are indicated by the blue boxes and stars. 



diative transition rates suggest the /i(1285), /i(1420) are 
somewhat flavor mixed as seen in the lattice data. 

Of particular interest to us here is the spectrum of ex- 
otic isoscalar mesons, where the basic structure of 1 ^ 
lightest below 0^ and 2^ observed for isovectors is 
again seen, but doubled to reflect the presence of both 
light and strange quarks. The positive parity exotics are 
very close to being ideally flavor mixed, but the 1 '", 
which in the isovector case we placed in a different super- 
multiplet, has a noticable degree of light-strange mixing. 
The large amount of strange-quark in the lighter state is 
most likely the dominant reason for it being significantly 
heavier than the corresponding isovector state^. 

In Figure [7] we show overlap histograms for the 
isoscalar 1 states. We see that the ^5*1 < ^5*1 < ^Di < 
^hyb^ pattern is again present, but doubled for light and 
strange quarks. The hybrid states show a greater degree 
of light-strange flavor mixing than do the qq states. Sim- 



* Note that suggestions from SU{3) Yang-Mills calculations are 
that the lightest exotic J^'^ glueball is much heavier than these 
states [33]. so we are justified in proceeding with a hybrid meson 
interpretation. 



ilarly we identify that the 2 + states near 2.4 GeV have 

[21 

strong overlap with the pnr x Dj^^ operators. Signals 
in the ^ sector are not clear enough to identify an an- 
ticipated pair of hybrid states above 2 GeV that would 
complete the hybrid supermultiplet. 

In summary we may well be observing the same hybrid 
supermultiplet structure in the isoscalar sector, and there 
appears to be some significant flavor mixing away from 
pure -y= (uu + dS) , ss structure. 



V. COMPARISON WITH MODELS 

The models described earlier can now be examined 
in light of the hybrid spectrum extracted from an in- 
terpretation of lattice QCD. We begin with the flux- 
tube model which predicts a large lightest supermulti- 
plet: (0,1,2)-+,1~,(0,1,2)+-,1++. The first four of 
these states have been clearly identified in the lattice 
QCD isovector spectrum, and one might argue that there 
are candidates for the remaining four at slightly higher 
energy. There would need to be considerable splitting 
within the supermultiplet and since two 2^ states ap- 
pear to be present in the mass region under consideration. 
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for light-quark overlap and dashed bars for strange. 



there would need to be a further flux-tube excitation al- 
most overlapping with the first. This appears to be an 
inefficient modeling of the observed spectrum. 

Bag models having a lowest energy "TE" gluon mode 
are chromomagnetic in gluonic character and give rise 
to a lightest hybrid supcrmultiplet that agrees with the 
one we observe. In |19| perturbative computations sug- 
gest that mixing of the non-exotic hybrid states with qq 
states should be small. There is some evidence in the 
lattice calculation (presented in Figure [4| for mixing for 
the ^ and 2 ^ states, but we have not attempted here 
to quantify the degree. Heavier hybrid meson states were 
not studied in detail in the bag model. 

p c 

Models in which a constituent gluon (Jg" " = 1 ) 
is added in S'-wave to a qq system completely fail to 
describe our observed states, giving rise to non-exotic 
positive parity states with quarks in an S'-wave and ex- 



otic and 1 with quarks in a P-wave. The 
is very high-lying in our extracted exotic spectrum and 
this suggests that this kind of hybrid construction does 
not contain the lowest energy gluonic excitation. 

If, one the other hand, we place a constituent gluon in 

_ PC 

a P-wave relative to the qq, and couple so that Jg" ° — 
1"' , we appear to be able to successfully describe both 
the lightest hybrid supcrmultiplet of (0, 1, 2) '", 1 (by 
having qq in an S'-wave) and the heavier exotic states, 
0"' ,(2"' )^ (with qq in a P-wave). In Appendix B we 
will present an explicit constituent gluon state construc- 
tion that can describe many of the overlap patterns ob- 
served. 

We should mention that another possible extension to 
the qq model of hadrons calls for the presence of states 
dominated by qqqq configurations and simple constituent 
quark mass counting suggests that they may enter at an 
energy scale similar to hybrid mesons. The main evidence 
against the importance of such constructions is the lack of 
any clear signals for exotic flavor mesons. Nevertheless it 
is important to include interpolators with good overlap 
onto such quark configurations into lattice calculations 
and appropriate constructions are currently underway for 
use in future computations. 



VI. HYBRID MESON PHENOMENOLOGY 

In this section we wish to address the current and pos- 
sible future experimental situation in light of what we 
have discovered about hybrid mesons in QCD. This will 
necessarily be somewhat speculative, as we have only de- 
termined the spectrum at unphysically heavy light quark 
masses, but we suspect that our qualitative observations 
will likely be robust under further reduction of TOtt. 

There is some evidence for states that might be candi- 
dates for our lightest isovector supcrmultiplet [3]. An ex- 
otic TTi near 1.6 GeV could be partnered with the claimed 
second excited 7r(1800) and a first excited 7r2(1880) to 
make up the full set of spin-triplet states. 

There is no clear candidate for the p-like member of 
the supcrmultiplet. This spin-singlet hybrid state would 
not be expected to be produced in e+e~ collisions, unless 
it mixes considerably with a qq^Si basis state. Figure |4] 
suggests it may not, but we must not infer too much from 
that plot. It presents the overlap for a particular smear- 
ing function that may accidentally suppress the overlap 
of a radially-excited state - some evidence for this kind of 
behavior is seen in Appendix A where a range of smear- 
ing functions are evaluated. Future lattice calculations 
should compute directly the vector decay constants of 
excited states {ipj^ip overlap with unsmeared fields, suit- 
ably renormalised) to more directly determine the possi- 
bility of production in e^e~. Other production processes 
such as photo- or electro-production off a hadron target 
may be more efficient than e'^e~ and again here one can 
in future use a lattice calculation to aid understanding. 
To the extent that photo- and electro- production can be 



11 



modeled by t-channel meson exchange, radiative transi- 
tion matrix elements, (M'| Vj7''-0|M), control the rate of 
production. These quantities can be extracted from lat- 
tice calculations, even for excited states, although to date 
only the charmonium sector has been considered .34] . 

Recently [35], a new isoscalar 1 state, F(2175), has 
been observed in e^e~ decaying into 0/o(98O) and (jyrj. 
This at first sight serves as a possible ss hybrid candi- 
date, but production in e+e" and the decay into (j) + X 
suggest quark spin-triplet structure (with the theoretical 
prejudice that quark spin tends to be preserved in such 
a decay). This is in disagreement with the quark spin- 
singlet structure of hybrid vector states that we have ob- 
served in this study. One way out of this is to note that 
the partial decay fractions io (f) -\- X have not been mea- 
sured and may in fact make up only a small contribution 
to the total width. The appropriate spin-singlet decays 
would be + X, but since almost nothing is known 
about the candidate state /ii(1380), observing these de- 
cays is not straightforward experimentally. 

One reason to propose hybrid character for a state is 
to argue that there is an overpopulation of states with 
respect to qq expectations, as has happened with the 
1^(4260) in the charmonium region. Such an argument 
does not hold water here as there is only one identified 
(f) excitation, the 0(1680), while we might expect at least 
two: a radially excited ss^Si and ss^Di. We may have 
to content ourselves with a more mundane explanation of 
the 1^(2175), certainly there is currently no overwhelming 
reason to assume it is a hybrid meson. 

In our lattice calculation of isoscalar mesons we found 
that the 1 members of the lightest hybrid supermul- 
tiplet may have a significant mixing of (uu + dcTj and 
ss, in contrast to the conventional vector states which 
appear to be very close to ideally flavor mixed. We need 
to verify that this property remains true as one reduces 
the quark mass, but it may eventually be useful as a 
phenomenological filter for hybrid character in the ex- 
perimental isoscalar vector meson spectrum. 

VII. SUMMARY 

We propose that, using the results of lattice QCD 
computations, we have identified the isovector mem- 
bers of lightest supermultiplet of hybrid mesons embed- 
ded within a spectrum of conventional mesons. The 



supermultiplet has quark spin-triplet states of J^^ = 
(0,1,2) and a quark spin-singlet state 1 . The 
identification followed from the presence of an exotic 
J^'-^' = 1 ^ state, the lack of convincing conventional 
qq supermultiplet (s) to describe the degeneracy pattern 
of the other states and the common overlap of these four 
states onto composite QCD operators having essentially 
gluonic structure. 

The form of the gluonic structure present in the oper- 
ators having good overlap with these states is chromo- 
magnetic, having Jg" " — 1^ . With the qq pair in an 
internal S'-wave this describes the observed J^'~^ . Heav- 
ier exotic hybrids 0"' , (2^ )^ are proposed to be qq P- 
waves coupled to the same gluonic excitation. A much 
heavier exotic state likely indicates the next distinct 
gluonic excitation form which might have Jg" " = 1 . 

The lattice calculations interpreted to come to the 
above conclusions are performed with light quark masses 
heavier than those realised in nature. Over the limited 
range considered (corresponding to pion masses between 
400 MeV and 700 MeV) the qualitative identification of 
the lightest hybrid multiplet appears to be robust. Fu- 
ture calculations will explore lighter quark masses. 

The nonexotic J^*^ hybrid mesons will not be eas- 
ily identified in experiment unless they happen to have 
some characteristic hadronic decays as suggested in cer- 
tain models [Sl|Tnil35'- In order for a lattice QCD calcula- 
tion to address the hadronic decays of hybrid mesons one 
requires there to be inclusion into the basis of operators 
some resembling pairs of mesons. Work in this direction 
is beginning with the mapping out of the amplitude for 
TTTT scattering in the non-resonant isospin-2 [37] channel 
and the resonant isospin-1 channel that contains the p 
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Appendix A: The effect of quark-field smearing in a bound-state model 



Within a constituent quark model, we propose a decomposition of the distillation-smeared fermion field as 



(27r)= 



^ 5,(|fc|) e-"=^u,(fc)bUfc) + e^'"i'.(fc)d^^(fc) 



(Al) 



where bl^(fc) annihilates a constituent quark of color i, spin a and momentum k and d||(fc) creates a constituent 
antiquark. The transformation that turns current quarks into constituent quarks will be not be specified here, we 
shall only assume that it is linear - an example would be the Bogoliubov transformation. We can allow there to be 
non-trivial energy-momentum dispersion, which might be expressed as a 'running' constituent quark mass, TOg(|fc|) 
featuring in the spinors. 

The standard implementation of distillation includes, with equal weight, a fixed number of eigenvectors of the 
gauge-covariant laplacian, — V^^„ — A„^„, as "Ylin ^n^ln which amounts to including all modes below some eigenvalue 
cutoff. We will include this in our model as a cutoff in quark momentum space, 5g(|fc|) = 0(A — jfcj). For a non- 
interacting theory in which the eigenvectors are plane waves ^'(a;) oc e^^'^ <5i,RGB subject to cubic boundary conditions, 
this modeling would be essentially exact. In the interacting theory there is still a similarity between the extracted 
eigenvalue distribution to the momenta allowed to a free colored particle in a cubic box as seen in Figure [8j 

Fermion bilinears without derivatives within this model will create constituent qq states from the vacuum 



d^x^'Ti^\x,t = 0)|0) = 



d^k 



E 



(2^)3 ^ J (2^)3 



d^k 



J2 S,i\k\)S,{\k\) u,{k)Tv^ik) d'xe'^^+^^y^ b'^{k)d'^{k)\0) 



/d^k ^ ..-^ ^ ^ ^ ^ 

j;^[S,{\k\)r Mk)rv,{-k) q'Ak)U-k) 



(A2) 
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FIG. 8: Grey points: Eigenvalues A„ of —V on a 24"^ dynamical lattice with ~ 400 MeV. Colored Boxes: k = (^) {ni + 
Tiy + n^) - allowed plane- wave momenta in a cube of side- length L — 24. 
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FIG. 9: Upper panel: Ratio of 1 first excited state to ground state overlaps for "/Onh" operator as a function of smearing 
radius, a, in spatial lattice units. Lower panel: Corresponding spectrum of lowest three 1~~ states. 



SO that rest-frame qq meson constructions of the following type will have overlap, 

\qq{n^'^+'Lj),mj)^ ^ (im^; 5ms| Jmj) ^ (is; is|5ms) / ^„l(|^|) M^) j^Ufc) (A3) 

Here (p is the radial momentum-space wavefunction. Which S", L, J combinations have non-zero overlap are determined 
by which T is used and can be worked out by expressing the spinors in the Pauli-Dirac basis. In all non-zero cases 

erlan nronortional to a, radial intepfral of the form 



one 



has the overlap proportional to a radial integral of the form 

k'dk [Sq{k)Yk'' <^„i(fc). 



(A4) 



and this integral can be accidentally small if Lp changes sign in the region of integration. We can explore whether 
this possibility manifests itself in practical calculations by varying the distillation smearing. A possible modified form 
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for the distillation smearing operator acting on quark fields is ^„ e '^^-^"/^ 6(A,„ax — A„) ^n^t where — V^^„ = A„^„ 
and where our default choice is cr = 0. In our bound-state model we will describe this via a smearing function 

= e-"''='/4e(A-A:). 

In a lattice QCD calculation we computed the spectrum of 1 isovector states with 700 MeV using a large 

operator basis for a range of smearing radii, cr, and the lowest 64 eigenvectors of — on a 16'^ lattice. In Figure|9]we 
show the low-lying spectrum extracted, which does not change qualitatively with cr, along with the ratio of overlaps 
onto pnr for the first-excited state and the ground state. We clearly see that for a particular region of smearing radii, 
cr 4 — 5, the overlap of the first excited state is suppressed onto an operator that accurately characterises its angular 
structure. This is an indication that we must be careful not to rule out a particular state assignment solely on the 
basis of small overlap values evaluated for one particular smearing function. 



Appendix B: Hybrid meson overlaps in a non-relativistic constituent gluon model 

. . _ . . [21 . 

We shall show using an explicit qqg state construction that our operators featuring Dj^^ can interpolate hybrid 



states featuring a J, 



1+ giuonic state. In order to consider the role of hybrids in the spectrum we must 



take account of the gauge-field part of the covariant derivatives. As mentioned in the text, the links that enter into 
discretised derivatives on the lattice are smeared. We will not assume that the corresponding smeared giuonic field 
is linear in the creation/annihilation of constituent gluons owing to its non-linear construction and the lack of an 
explicit transformation from Lagrangian gluons to constituent gluons. Our assumed form will be 



<Pq V 

(2^)3 



SM) 



©(a^a^) +0(aa) 



We will not have cause to go beyond one power of the constituent gluon creation operator in our interpretation 
of the spectrum. The constituent gluons created by a^^(g) have an energy-momentum dispersion that is not 

specified, but it is assumed to feature a mass-gap as |(7| — )• 0. Since stout-link smearing [3T] doesn't change the gauge- 
transformation properties of the links we'll make the simplifying assumption that the smearing function, 5g(|/c|) is 
color- independent and simply acts to damp out high-momentum modes. We have assumed that the constituent gluons 
have only transverse helicity, A = ±1, as would be appropriate in Coloumb gauge. The only piece of this decomposition 
that we will make use of in interpolating hybrid meson states is the term linear in the creation of a constituent gluon. 

We shall work under the assumption that a mass gap in the constituent gluon dispersion is such that states with 
more than one constituent gluon appear high in the spectrum and can be neglected. We construct qqg states at rest 
non-relativistically as in |21j but couple the internal angular momentum in a slightly different order, 



qq{n^'^'LjJg[j^^]- J 



PC 



mj 



(27r)3 (27r)3 ^Vi'^i'i.-i; ^ 
X (im; ^m\Sms){LmL; Sms\Jqqmqq){jqqmqq; Jgmg\jmj) 



X 2'-|j 



(Bl) 



where the sum is over all z-components of angular momenta and helicities, rn,fn,rnL,ms,mqq,rng, fi and where 
^ = Pg{—l)^'s. The form of the dependence on the magnitudes of quark and gluon momenta, "if, will not be specified 
here. The presence of a generator of SU{3) color, t?j ensures the qq pair are coupled to a color octet. A relativistically 
invariant construction would need to account for Wigner-rotations arising from the qq not being at rest within a qqg 
meson and while these constructions can be performed, |39) . they are somewhat elaborate and we will not attempt 
them here. 

Operators with a single derivative (-Dj'^j^) are capable of interpolating qqg states through the A piece of the 



derivative, D 



— > 



d — 2igA as well as qq states through the conventional derivative. For qqg the overlap will be 



[qqln^'^+^LjJglJ^^]; J— ,m.j 



tPC 



r (fp (fq 
°^ J (27r)3 (27r)3 ^ ' 
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Jd = 1 


Jd = 2 
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Ifcpy--(fc) 


[DA]{qqg) 
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(Imi; lm2 2mD)Pmiem2('?i'^) 


[F](gw) 





(gxr(g,A))^^ 






TABLE II: Angular integrand in overlap of qq and qqg with D 



[2] 



constructions (as in 



B3l . 



To be somewhat consistent with our neglect of Wigner-rotations in the qqg state constructions we will not consider 

terms in the spinor bilinear uTv that are proportional to the total qq momentum, —q. In this case the integral over q 

can be performed and is non-zero only if Jg = 1 and rUg = m^). The resulting overlap is proportional to 1 + ^ = 1 — Pg 

p 

so that the gluonic excitation interpolated has Jg" = 1~. What remains in the overlap are radial integrals and the 
integral over the qq relative momentum direction, p, which determines which qq L, S can contribute. 
With two derivatives we can consider four possibilities. 



— > 



- 2ig{ dmiA 

r. 

+ 2ig{d,n^A 

rr. _ 



Ar 



A 

m2 -^m-i 



i5 [Am, Alls]) 



[DD] 
[DA] 
[F] 
[AA]. 



(B3) 



In Table ITT] we present these possibilities projected into definite a ngu lar momentum for qq and qqg. Performing the 

— [21 . P - 

we find that D'J^^Q^^ can interpolate Jg " = 1 

[21 P 

gluonic states while Dy^^^ interpolates Jg ' 



Bl 



angular integrals in an overlap with states constructed as in Eqn. 

1+. The Clebsch-Gordan couplings are such in this model that the 

[21 

overlaps for the lightest hybrid supermultiplet onto the operators {pnRi '"'nr} ^ jLi would be expected to have the 
same value 



'S,)g[l 



+1 



0" 



qq{^Si) g[l+]; r 
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(B4) 



As seen in figure [3] this is observed approximately in the lattice data for the states suggested to form the lightest 
hybrid supermultiplet. 
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